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SECTION  I 

INTRODUCTION  AND  SUMMARY 


In  recent  years,  the  Rotor-Bearing  Dynamics  Design  Technology  Series, 
AFAPL-TR-65-45  (Parts  I through  X)  has  been  considered  by  many  engineers 
to  be  an  important  part  of  their  basic  analytical  tool  kit.  However, 
since  the  issuance  of  the  first  volume,  in  May  of  1965,  the  state-of-the- 
art  has  significantly  advanced.  Further,  new  techniques  of  data  presen- 
tation have  been  developed,  computer  capabilities  have  increased  and  some 
minor  typographical  and  technical  errors  were  uncovered. 

/Part  IV  of  AFAPL-TR-65-45  treated  design  data  for  typical  deep-groove 
and  angular  contact  bearings.  The  data  was  presented  in  graphical  form 
and  consisted  of  direct  radial  stiffness,  load  carrying  capacity,  and 
load  levels.  In  addition  design  guidelines  and  limitations  were  discussed. 
The  major  deficiencies  of  this  original  volume  were  that  centrifugal  effects 
due  to  high  speed  were  ignored,  and  axial  and  angular  stiffness  informa- 
tion were  omitted. 

Subsequent  to  the  publication  of  Part  IV,  several  extensive  treatments  of 
ball  bearings  including  el astohydrodynamic , thermal,  and  cage  effects  have 
been  published.  The  computer  program  of  Mauriello,  LaGasse,  and  Jones  (3) 
considers  both  elastohydrodynamic  and  cage  effects.  The  more  recent  computer 
based  design  guide  prepared  by  Creceliuc  and  Pirvlcs  (4)  treats  elastohydro- 
dynamic, thermal,  and  cage  effects  for  a system  of  ball  and  roller  bearings. 

Thus  very  sophisticated  analytical  tools  are  available  for  the  design  and 
application  of  ball  bearings.  Neither  of  these  tools,  however,  provide 
the  user  with  the  stiffness  matrix  required  for  solution  of  rotor  dynamics 
problems.  In  addition  both  computer  programs  are  very  large  and  require  an 
extensive  computer  facility  for  use. 

The  present  volume  is  intended  as  an  update  of  the  original  Part  1V(1). 

Those  aspects  of  the  original  Part  I V ( 1 ) which  treated  general  design 
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aspects  of  ball  bearings,  load  capacity,  speed  limitations,  etc.  have  been 
deleted  since  their  coverage  is  superficial  compared  to  the  more  sophisti- 
cated computer  tools  now  available  (3,4).  Only  those  parts  directly  con- 
nected with  preparation  of  input  for  the  rotordynamic  response  programs 
(Volume  I of  the  revised  series)  have  been  retained.  The  complete  stiffness 
matrix  is  calculated  including  centrifugal  effects.  Considerations  such  as 
elastohydrodynamic  and  cage  effects  are  not  included  since  they  have  little 
influence  on  the  calculation  of  ball  bearing  stiffness.  The  resulting 
program  (Appendix  A)  is  reasonably  small  and  easy  to  use.  Lastly,  the  stiff- 
ness data  included  in  the  original  Part  IV  have  bean  updated  and  are  included 
in  Appendices  B,  C and  D. 


SECTION  II 


ANALYSIS 


2.1  General  Bearing  Model  and  Coordinate  System 

Accurate  calculation  of  the  lateral  dynamic  response  of  a high-speed  rotor 
depends  on  realistic  characterization  of  the  support  bearings.  In  the  most 
general  case,  both  linear  and  angular  motions  are  restrained  by  the  support 
bearings  at  the  attachment  location.  In  the  analytical  model,  the  reaction 
force  and  the  reaction  moment  of  each  bearing  are  felt  by  the  rotor  through 
a single  station  of  the  rotor  axis.  As  schematically  illustrated  in  Figure 
la,  a coil  spring  restraining  the  lateral  displacement  and  a torsion  spring 
which  tends  to  oppose  an  inclination  are  attached  to  the  same  point  of  the 
rotor  axis.  A complete  description  of  the  chr.racteristics  of  the  support 
bearings,  however,  involves  much  more  than  the  specification  of  the  two 
spring  constants.  This  is  because: 

• The  lateral  motion  of  the  rotor  axis  is  concerned  with  two 
displacement  components  and  two  inclination  components. 

• The  restraining  characteristics  may  include  cross  coupling 
among  various  displacement/inclination  coordinates. 

• The  restraining  force/moment  may  not  be  temporally  in  phase 
with  the  displacement/inclination. 

• The  restraining  characteristics  of  the  bearing  may  be 
dependent  on  either  the  rotor  speed  or  the  frequency  of 
vibration,  or  both. 

• Bearing  pedestal  compliance  may  not  be  negligible. 


To  accommodate  the  above  considerations,  the  support  bearing  characteristics 
are  described  in  Reference  2 by  a four-degrees-of- freedom  impedance  matrix  as 
defined  in  Equation  (1) : 
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(1) 


where  WV1  is  a 
— N 

displacements 


column  vector  containing  elements  which  are  the  two  lateral 

(6,6)  and  the  two  lateral  inclinations  (0  , 0 ) of  the 
x y x y 
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(a)  Bearing  Stiffness  Model 


(b)  Bearing  Location 
Coordinate  System 


Figure  1 
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rotor  axis  at  the  bearing  station  N. 

Employing  a right  handed  Cartesian  representation  in  a lateral  plane  as 

depicted  in  Figure  lb,  the  2-axis  is  coincident  with  the  spin  vector  of  the 

rotor.  The  x-axis  is  oriented  in  Che  direction  of  the  external  static  load, 

and  the  y-axis  is  perpendicular  to  both  z and  x axes  forming  the  right 

handed  triad  (x,  y,  z).  (6,6)  are  respectively  lateral  lineal  displace- 

x y 

mont  components  of  the  rotor  axis  along  the  (x,  y)  directions.  (8,9) 

x y 

are  lateral  inclination  components  respectively  in  the  (z-x,  z-y)  planes. 

Note  that  9 is  a rotation  about  the  y-axis,  while  9 is  a rotation  about 

x y 

the  negative  x-axis. 


is  a complex  (4x4  matrix),  and  in  accordance  with  the  common  notation 
for  stiffness  and  damping  coefficients,  may  be  expressed  as 

SN  * ^ + ^§N  (2) 


where  is  the  stiffness  matrix  and  is  the  damping  matrix.  v is  the'* 
frequency  of  vibration,  host  commonly,  lateral  linear  and  angular  displace- 
ments do  not  interact  with  each  other  so  that  the  non-vanishing  portions  of 
and  are  separate  2x2  matrices.  That  is 
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Accordingly,  a total  characterization  of  a support  bearing  would  include 
sixteen  coefficients  which  make  up  the  4 (2x2)  matrices: 
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(8) 


In  the  event  that  the  pedestal  compliance  is  significant,  then  the  effective 
support  impedance  can  be  calculated  from 

- 'S»  ■ ip'1-'  (9) 

where  subscripts  "p"  and  "b"  refer  to  the  pedestal  and  bearing  respectively. 
Note  that  both  pedestal  inertia  and  damping  may  be  included  in  2^ . 

2.2  General  Bearing  Support  Characteristics 

The  function  of  a bearing  is  to  restrict  the  rotor  axis  to  a nominal  axis 
under  realistic  static  and  dynamic  load  environments.  Deviation  of  any 
particular  point  of  the  rotor  axis  from  the  nominal  line  can  be  character- 
ized by  three  lineal  and  two  angular  displacements.  These  may  be  designa- 
ted as  (6  , 6 , 6 , 9 , 6 ) in  accordance  with  a right-handed  Cartesian 
x y z x y 

reference  system.  The  z-coordinate  is  coincident  with  the  reference  axis 
and  is  directed  toward  the  spin  vector.  (0^,  9^)  are  rotor  axis  inclina- 
tions respectively  in  the  z-x  and  z-v  planes.  The  x-coordinate  is  directed 
toward  the  predominant  static  load;  e.g.,  earth  gravity.  Ideally,  the 
bearing  would  resist  the  occurrence  of  any  displacement  so  that  the 
reaction  force  system  imparted  by  the  bearing  to  the  rotor  is  generally 
expressed  in  matrix  notation  as 

F = -Z-x  (10) 
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F is  a column  vector  comprising  the  five  reaction  components  (F^,  F , F^, 

M , M ),  while  x is  the  displacement  vector  (6, 6, 6, 9, 9).  Zisa 
xy  — x y z x y 55 

(5  x 5)  matrix  containing  the  elements  with  both  indicies  (i,  j) 
ranging  from  1 to  5.  The  values  of  characterize  how  rotor  displace- 
ments are  being  resisted  by  the  bearing. 


From  the  standpoint  of  dynamic  perturbation,  distinction  is  made  between 
a static  equilibrium  component  and  a dynamic  perturbation  component  for 
both  the  displacements  and  the  reactions.  Thus, 

x - x^  +X'i  F = Fo+  F*  (11) 

(x1 . F' ) are  respectively  presumed  to  be  infinitesimal  in  comparison  with 

(x  , F ).  Accordingly,  Z..  are  regarded  as  dependent  on  x but  not  on  x'- 
-o  — o ij  — o ~ 

To  illustrate  the  idea  of  perturbation  linearization,  one  may  examine  the 
one-dimensional  load-displacement  curve  shown  in  Figure  2. 


As  illustrated,  the  load-displacement  relationship  is  a 3/2  power  law  in 
accordance  with  the  Hertzian  point  contact  formula.  It  is  not  possible 
to  describe  the  entire  range  by  a linear  approximation.  However,  if  a 
small  dynamic  perturbation  is  taken  around  a static  equilibrium  point, 

6'  < 6^  , the  small  segment  of  the  load-displacement  curve  can  be 

approximated  by  a local  tangent  line.  The  corresponding  force  increment 
is 


F' 

x 


(12) 


where  6'  is  the  incremental  displacement.  3F  /96x  will  depend  on  the 

amplitude  of  6^  . 

o 

The  question  of  history  dependence  is  resolved  by  regarding  x'  as  periodic 
motions  at  any  frequency  v of  interest,  and  Z accordingly  would  have 
both  real  and  imaginary  parts  and  may  also  be  depended  on  both  the  rotor 
speed  uj  and  the  vibration  frequency  v. 


To  avoid  notational  clumsiness,  the  primes  will  he  dropped  from  (F',x') 
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! m 


which  are  understood  to  be  dynamic  perturbation  quantities  unless  the  sub- 
script "0"  is  used  to  designate  the  static  equilibrium  condition. 


2.3  Ball  Bearing  Characterization 

In  many  ways  the  ball  bearing  is  much  simpler  to  model  from  a rotor  dynamic 
point  of  view  than  a fluid  film  bearing.  In  general,  the  followirg  two 
simplifications  can  be  made: 

• The  restraining  characteristics  do  not  include  cross 
coupling  among  the  various  displacement/inclination 
coordinates . 

The  restraining  force/moment  is  normally  temporally 
in  phase  with  the  displacement/inclinufion. 

Figure  3 shows  a ball  bearing  referred  to  in  an  orthogonal  xyz  coordinate 
system.  The  outer  ring  is  fixed  but  the  inner  ring  may  move  with  respect 

to  the  coordinate  system.  Both  rings  are  free  to  rotate  about  their  axes. 

.....  «♦  • 

Three  lineal  displacements,  , -5^,  and  two  angular  displacements , 0^, 

9 , are  required  to  define  the  spatial  position  and  attitude  of  the  inner 

V 

ring  when  it  is  displaced  from  its  initial  position.  For  purposes  of 
derivation  the  initial  situation  is  that  existing  when  the  bearing'?  end 
play  is  just  taken  up  in  the  thrust  direction.  Figure  3 shows  these 
displacements  in  the  positive  sense. 

Figure  4 shows  some  important  dimensions  and  establishes  the  convention 
of  the  ball-position  index  q.  The  contact  angle  f?  is  the  initial  mounted 
contact  angle  and  is  shown  in  the  positive  sense. 

2.3.1  Stiffness 

The  total  characterization  of  a ball  bearing’s  stiffness  can  be 
expressed  Kv  the  natrix  at  the  top  of  page  12. 
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(13) 


The  lineal  and  angular  stiffness  matrices  (Equations  5 and  7)  ran  be  derived 
from  Equation  (13).  For  example: 


(K) 


lineal 


(K> 


angular 
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(14) 


(15) 


Note  that  although  the  axial  components  of  stiffness  are  not  utilized 
by  the  lateral  rotor  dynamics  program  (2),  they  have  been  retained 
in  the  general  ball  bearing  stiffness  matrix.  Equation  (13).  The 
axial  stiffness  would  be  required,  for  example,  if  the  reader  was 
calculating  the  axial  natural  frequency  of  a b3ll  bearing  mounted  shaft. 


2.3.2  Damping 

There  is  very  little  data  on  ball  bearing  damping.  The  data  included 
in  Reference  1,  which  is  for  nonrotating  grease  packed  bearings, 
suggested  a value  in  the  order  of  15-20  pound  sec/in.  This  should  be 
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used  only  as  an  approximate  figure  since  this  has  not  been 
confirmed . 


2.4  Load-Deflection  Relationships 

For  a given  ball-raceway  contact,  the  load  deflection  relationship  is  given 
by  an  equation  of  the  form: 

3/2 

P - K AJ/  (16) 


The  total  normal  approach  between  two  raceways  under  load  separated  by  a 
rolling  element  is  the  sum  of  the  approaches  between  the  rolling  element 
and  each  raceway.  Hence, 


(17) 

(18) 

where  and  K2  are  a function  of  the  ball-race  geometry  and  material 
properties . 


A1  +A2 


and 


-.3/2 


(1/K| ) 2/3  -f  (1/K2)2/3 


2.5  Bali  Bearings  Under  Radial  Load 

For  a rigidly  supported  bearing  subjected  to  radial  load,  the  radial 
deformation  at  any  rolling  element  angular  position  is  given  by 

A - « cos*  - (1/2 )P  (19) 

* max  u 

in  which  6 is  the  maximum  deflection,  occurring  at  4 « 0 and  P_  is  the 
max  o-r  d 

diametral  clearance.  Figure  5 illustrates  a radial  bearing  with  clearance. 

Equation  (19)  may  be  rearranged  in  terms  of  maximum  deformation  as  follows: 

A " A U - ~ (1  - cos*)}  (20) 

max  l £ 

in  which 
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(21) 


It  is  clear  from  Equation  (21)  that  the  angular  extent  of  the  load  zone  is 
determined  by  the  diametral  clearance  such  that 


(22) 


For  zero  clearance,  $ * 90 
From  Equation  (16) 


- - 

— — 
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(23) 


Therefore,  from  Equations  (20)  and  (23) 


P 

max 


1 - (1  - cos*) 


3/2 


(24) 


For  static  equilibrium  to  exist,  the  applied  radial  load  must  equal  the  sum 
of  the  horizontal  (parallel  to  the  load  direction)  components  of  the  rolling 
element  loads: 
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(25) 
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1 - 77  (1  - cos<p* 


3/2 


cos<)> 


(26) 


The  summation  in  above  equations  applies  only  to  the  angular  extent  of  the 
load  zone.  Equation  (26)  can  also  be  written  in  integral  form: 
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or 
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in  which 
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(23) 


(29) 


The  radial  integral  of  Equation  (29)  has  been  evaluated  numerically  by 
Harris  (5)  for  various  values  of 


From  Equation  (16) , 


K A 


3/2  _ 


i l'  U , A 

max  <t>=0 


K { 6 - \ PJ3/2 

max  2 D 


Therefore, 


i 3/2 

F = nK  {6  - tt  Pn}  J (?) 

x max  2D  x 


(30) 


(31) 


For  a given  bearing  with  a given  clearance  under  a given  load,  Equation  (31) 
may  be  solved  by  trial  and  error. 

For  ball  bearings  under  pure  radial  load  and  zero  clearance,  ” .2285, 
and  it  can  be  shown  from  Reference  6 that 

4.37F 

P * -x-  (32) 

max  ncosB 


Accounting  for  nominal  diametral  clearance  in  the  bearing,  one  may  use  the 
following  approximation, 

5F 

P = S-r  (33) 

max  ncosB 


2.5.1  Approximation  to  Radial  Stiffness 

A complete  definition  of  bearing  stiffness  requires  calculation  of  the 
stiffness  matrix  defined  by  Equation  (13).  However,  for  many 
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applications,  simplifications  can  be  made  which  permit  calculation 
of  radial  stiffness  of  a ball  bearing  under  radial  load  alone. 


Palmgren  (6)  gives  a series  of  formulas  to  calculate  bearing  deflec- 
tion for  specific  conditions  of  loading.  For  slow  and  moderate 
speed  deep-groove  and  angular-contact  ball  bearings  subjected  to 
radial  load  which  causes  only  radial  deflection,  that  is  6 =0, 

p 2/3 

sx  ' i'58  * 10'5 -T77- - <w> 

d cosB 


For  self-aligning  ball  bearings, 

2/3 


6 = 2.53  x 10 

x 


, P 
-5  max 


a3-/3  a 
d costs 


(35) 


For  this  case  the  maximum  rolling  element  load  has  previously  been 
shown  to  be 


P 

max 


5F 

x 

ncosS 


(33) 


Substituting  for  Pmax  in  Equacion  (34)  yields 

4.62  x 10~5F  2/3 
* x 

x = 2/3.1/3  5/3 

n d cos  8 


(36) 


Transposing,  then  taking  the  derivative  with  respect  to  x,  the  stiffness 
of  the  bearing  is 

„ , ,,  in6  .1/2  5/2..  1/2 

K * 4.77  x 10  nd  cos  66  (37) 

XX  X 

For  self-aligning  ball  bearings 

V O .,1/2  5/2  1/2 

K * 2.36  x 10  nd  cos  66  (38) 

XX  X 

From  Equations  (37)  and  (33)  it  is  apparent  that  the  deflection-stiff- 
ness relationship  is  nonlinear  because  it  is  dependent  on  the  square 
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root  of  radial  deflection.  In  this  respect,  a ball  bearing  is 
unlike  a simple  spring  for  which  deflection  is  linear  with  respect 
to  load. 


By  substituting  for  6^  in  Equations  (37)  and  (38),  bearing  stiffness 
can  be  expressed  in  terms  of  load  as: 


Angular  Contact  Bearings 


xx 


.1/3  2/3  5/3. 

d n cos  8 


3.247  x 104  F 1/3 
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Self-Aligning  Bearings 


xx 


.1/3  2/3  5/3 

d n cos  8 


4 1/3 

2.028  x 10  F 
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These  expressions  are  plotted  in  Figure  6. 


(39) 


(40) 


2.6  Bearings  Under  Thrust  Load 

Thrust  ball  bearings  subjected  to  a centric  thrust  load  have  the  load 
distributed  equally  among  the  rolling  elements.  Hence, 


P 


F 

z 

n sin8 


(41) 


In  Equation  (41),  8 is  the  contact  angle  which  occurs  in  the  loaded  bearing. 
For  thrust  ball  bearings  whose  contact  angles  are  nominally  less  than  90 
degrees,  the  contact  angle  in  the  loaded  bearing  is  greater  than  the  initial 
contact  angle  8'  which  occurs  in  the  unloaded  bearing. 


In  the  absence  of  centrifugal  loading,  the  contact  angles  at  inner  and 
outer  raceways  are  identical;  however,  they  are  greater  than  those  in  the 
unloaded  condition.  In  the  unloaded  condition,  contact  angle  is  defined  by 


cos6' 


(42) 
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in  which  is  the  mounted  diametral  clearance.  A thrust  load,  Fz , applied 
to  the  inner  ring  as  shown  by  Figure  7 causes  an  axial  deflection  6 . This 
axial  deflection  is  a component  of  a normal  deflection  along  the  line  of 
contact  such  that  from  Figure  7 


A 


Bd  ( 


cosB 1 

cosB 


1) 


3/2 

Since  P * KA  , 
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K(Bd)3/2 


cosg ' 
cosB 


(43) 


(44) 


Substitution  of  Equation  (41)  into  Equation  (44)  yields 


nK(Bd) 


3/2 


cosB'  3/2 

sln6  Cr  - » 


(45) 


Since  K is  a function  of  the  final  contact  angle,  8,  Equation  (45)  must  be 

solved  by  trial  and  error  to  yield  an  exact  solution  for  8.  Jones  (7)  has 

defined  an  axial  deflection  constant  K as  follows: 

z 


K 

z 


B 

g(-hr)  + g(-y) 


(46) 


in  which  y 3 - — °S  ^ and  g(+y)  refers  to  the  inner  raceway  and  g(-y)  refers 

to  the  outer  raceway.  Jones  further  indicates  that  the  sum  of  g(+y)  and 

g(-y)  remains  virtually  constant  for  all  contact  angles  being  dependent  only 

on  total  curvature  B.  The  axial  deflection  constant  K is  related  to  K as 

z 

follows : 


K i112 

1 57T 


(47) 


Hence : 

z 


sing 


,cosB 1 
cosB 


(48) 


Taking  K from  Reference  7,  Equation  (48)  may  be  solved  numerically  by  the 
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Newton-Raphson  method. 


The  axial  deflection  6^  corresponding  to  A may  also  be  determined  from  Figure 
7 as  follows: 


(Bd  + A)  sing  - Bd  sind’ 


Substituting  A from  Equation  (A3)  yields 


Bd  sin  (6-3’) 
cosd 


2.6.1  Approximation  to  Axial  Stiffness 


From  examination  of  Figure  7,  the  following  approximation  can  be  made: 


A ~ 6 sir 8 
2 


Substituting  for  A in  Equation  (16)  gives 


1 P‘ 


K273  sing 


Using  the  Jones  (7)  axial  deflection  parameter,  Equation  (32)  becomes 

,,2/3 


Z K 2/3d1/3  s in(3 
2 

A similar  approximation  was  arrived  at  by  Palmgren  (6).  Suggested 
values  for  axial  deflection  under  pure  axial  load  are: 
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6 = 2.53  x 10 

2 


.1/3  . a 
u sing 


22 


Thrust  Ball  Bearings 


£ 

z 


1.9  x 10 


-5 


?2./3 

' max 


d1/3sine 


The  corresponding  stiffnesses  are: 


Angular-Contact  Bearings 
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zz 
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9.5  x 104  F1/3 
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(56) 


(57) 


Self-Aligning  Ball  Bearings 
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.1/3  2/3  , 5/3. 
d n sin  6 


= 5.9  x 104  F1/3 
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(58) 


Thrust  Ball  Bearings 
K 


zz 


4 a/3 


d1/3n2/3sl„5/3S 


= 7.9  x 10  F 


(59) 


These  relationships  are  plotted  in  Figure  8 cor  the  three  types  of 
bearings . 

2.7  Ball  Bearings  Under  Combined  Loading 

Except  for  ball  bearings  under  simple  radial  or  thrust  load,  treated 
previously,  there  are  very  few  solutions  for  bearing  stiffness  that  can 
be  evaluated  by  simple  hand  computation. 

When  a ball  bearing  operates  at  high  speed,  the  body  forces  resulting  from 
the.  ball's  motion  become  significant  and  must  be  considered  in  any  analysis. 
Figure  9 shows  the  forces  and  moments  acting  on  the  ball  in  a high-speed 

ball  bearing.  The  operating  contact  angle  at  the  outer  contact  is  less  than 
that  on  the  inner  because  of  the  body  forces.  Subscript  1 refers  to  an 
outer  contact  and  subscript  2 to  an  inner. 
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Thrust  Stiffness  Parameter  vs.  Thrust  Load 


1 


The  ball 


In  Figure  10  the  ball  center  is  fixed  in  the  plane  of  the  paper 
rotates  with  the  angular  velocity  u directed  at  the  angle  a.  . race 

n 

rotates  about  the  bearing  axis  with  the  angular  velocity  w , i ,2  relative 
to  the  retainer.  For  the  linear  velocity  of  race  and  ball  to  be  equal  at  a 
contact,  the  following  relation  must  be  satisfied: 

C , £'  (1  + C -y  )u> 

“B  = dco~s(S^-a)  <60) 

where  = 1 and  C2  = -1 
dcosS^ 

Yi  " ~Er~~  (61) 

For  the  iC^  race  to  be  stationary,  the  ball  must  orbit  with  the  angular 
velocity  such  that 

aE  = ""i  (62) 

With  stationary  outer  race  and  rotating  inner,  the  actual  angular  velocity  of  the 
inner  race  is 

^2  ” >J2  + = "°2  ~ Jl  (63) 

For  rotating  outer  and  stationary  inner,  the  actual  angular  velocity  of  the 
outer  is 


*”  *^1  ^ " W ^ o>2 

From  Equation  (60), 

(1-^2)003(6^0) 

^2  (l+y^)cos  (62-01) 

considering  that  both  races  may  rotate  there  results: 

(^1-^2 ) (1-Y2^cos  (3j_~a) 
jJ1  (l+Y1)co8(62*-n)  + (l-y2)cos(B1-a) 


(64) 


(65) 


(66) 
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Figure  10  Ball  Motion  Vectors 
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(67) 


(l+Y^cos^-ct) 

w'2  (l+Y1)cos(62-a)+(l-y2)cos(61-a) 

E'{(n1-n2)(i+Yl)d-Y2)} 

UB  d{  (l+Y1)cos(S2-a)  + (l-Y2)co8(B1-c») } 

(1+Y1)cos(62~^)+‘72  (l-Y2)C0S(6i-a) 
E (l+Y1)cos(B2-a)+(l-Y2)cos(81-a) 


(68) 


(69) 


For  an  arbitrary  choice  of  a there  will  be  a spin  of  the  ball  relative  to  a 
race  about  the  normal  at  the  center  of  the  contact  area.  This  is  illustra- 
ted in  Figure  11. 

From  Figure  10, 

u>  ■ C.{-a).sinS.  + u)nsin(6  -a)}  (70) 

i i i B i 

The  controlling  race  hypothesis  assumes  that  all  spin  occurs  at  one  contact 
while  no  spin  occurs  at  the  other.  The  contact  at  which  no  spin  exists  is 
called  the  controlling  race.  Lightly-loaded  bearings  may  depart  somewhat 
from  this  situation. 


If  u)a^  is  made  zero  and  Equation  (70)  solved  for  a,  there  results  for 
outer  race  control: 

sinB.cosB 

a * tan j (71) 

cos  6^+v  ^ 

and  for  inner  race  control: 

sinBpCosBp 

a = tan-1  j- (72) 

cos  $2~'*2 

The  existence  of  a particular  type  of  control  depends  on  the  relative  torques 
required  to  produce  spin  at  the  two  contacts. 


The  torque  required  to  produce  spin  is: 


e 


Figure  11  Illustration  of  Race  Control  Concept 


0 = 3uP  a E(e) 

4 8 

Q = spin  torque,  lb-in 

P = contact  load,  lbs. 

u ■ coefficient  of  sliding  friction 

a * semi-major  axis  of  pressure  ellipse,  in. 

E(e)  = complete  elliptic  integral  of  the  second  kind 

formed  with  the  modulus  sine  where  cose  = b/a, 
b being  the  semi-minor  axis  of  the  ellipse 


(73) 


Figure  12  shows  a ball  acted  on  by  the  spin  torque  vectors  and  Outer 

race  control  will  exist  if  projected  on  0^  is  greater  than  c^.  Inner 

race  control  will  exist  if  projected  on  is  greater  than  . 


Figure  13  shows  the  relative  positions  of  ball  and  race  curvature  centers 
before  and  after  the  application  of  the  five  displacements,  5^,  6^,  6^, 


9 , and  9 . 

x y 


The  outer  race  curvature  center  is  fixed.  The  inner  race 


curvature  center  has  moved  to  (A^.A^)  and  the  ball  center  to  (X^^). 

A - Bdsin?  + A + A"  + r[  (9  + 9")  8in4>  + (9  + 9")cosd>]  (74) 

1 7.  7.  X X V V 


A?  *>  BdcosB  + (6  + 5")  cos}  + (6  + 6")  sin$  — — 

z xx  y y z 


(73) 


B - f + f - 1 


R = E/2  + ( f 2 - -5)  dcosB 


(76) 

(77) 


PD 

s 

diametral  clearance, 

in . 

6" 

X 

B 

initial 

displacement 

along 

x, 

in . 

6:| 

y 

- 

initial 

displacement 

along 

y> 

in. 

6" 

z 

B 

Initial 

displacement 

along 

z, 

in. 

9” 

X 

S 

initial 

misalignment 

about 

X, 

rad . 

9" 

V 

S 

initial 

misal ignment 

about 

y, 

rad . 

For  the  forces  and  moments  acting  on  the  ball  to  be  in  equilibrium,  the 
following  must  be  satisfied: 


P^sinB^ 


2M 

P-sinB-  - ~ { XcosB. -(l-\)cos8.) 
/ la  1 L 


(78) 
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Figure  12  Spin  Torque  Vectors 
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FINAL  POSITION  OF  INNER 
RACE  CURVATURE  CENTER 


Figure  13 


Race  Curvature  Center  Deflection 
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PjCOsBi  " P2co8®2  + "d”  {AsinB^-d-X^inBj)  - F£  * 0 


where  X is  1 for  outer  race  control  and  0 for  inner  race  control  and 

F ■ M ftp  e'/2 
c t 


l?  nE  WB  8ina 


E'  » Operating  Pitch  Diameter,  in. 
- E + 2X2  - (f  - . 5)dcos3 


Prom  Figure  13 


, 'I  ^ 

tan  — 


^2 

-1  Vi 
A2"X2 


A^  and  A2  are  the  elastic  approaches  of  ball  and  raceways  at  the  contacts 

A1  “ (\2  + X22)1/?  ~ (fi  ~ -5>d  > 0 (85) 

2 2^/2 

A2  - {(Ax  - Xj_)z  + (A2  - X2)^}  - (f2  - . 5)d  > 0 (36) 

The  contact  loads  are  related  to  the  elastic  approaches  through 

Pi  * Ki  Ai/2  (37) 


where 


11.34771 


E(ci)d 


(np  + nR)cose  3 

B 1 {K(et)}  (4  - l/fi 


4 (1~VR  > 


2Vi 

1+Vi 
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(90) 


ui-V) 


E_  and  E_  are  moduli  of  elasticity  for  race  and  ball. 

K 5 

v_  and  v_  arc  Poisson's  Ratios  for  race  and  ball. 

K D 

K(e)  and  E(e)  are  the  complete  elliptic  integrals  of  the  first  and  second 
kind  having  the  modulus  sine. 


sine 


'1  - 


(91) 


b^  and  a^  are  the  semi-minor  and  setai-major  axes  of  the  pressure  ellipse 
e is  related  to  the  contact  geometry  through 


9 (1-cost . )E(e . ) 

cot“e1  = ,_ 


2{K(ei)  - E(e.)} 


(92) 


COST 


2Vi 

1+CiYi 


1 2C.y« 

4 _ J_  _ _ L-L. 

fi  1+CiYi 


(93) 


Equations  (73)  and  (79)  are  a set  of  nonlinear,  simultaneous  equations.  When 
the  relative  ring  displacements  are  held  constant,  the  variables  in  Equations 
(78)  and  (79)  are  and  These  are  evaluated  numerically. 

The  reactions  of  the  bearing  on  the  shaft  are 


F' 

x 


n 2 ( 1—  X )M 

I [P2  cose  + — J- — 1 sine2  3 cos* 
1=1  i i i 


(94) 


F' 

y 


n 2( l-X ,)M 

n i g< 

£ [p~  cos8?  + ~ 

i=l  i i 


j sinS0  ] sin^ 


(93) 


F' 

z 


2 ( l-X . ) M 

n i.  g . 

I [P-  sin8~  — cosB9  ] 

i=l  i i d “i 


(96) 
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fjdjJsin^  (97) 

f2<l)]cos^i  (98) 

Then  the  equilibrium  of  external  forces 

l {F*  + F } = 0 (99) 

i»x,y,z 

F^  - external  load  along  x,  lbs. 

Fy  - external  load  along  y,  lbs. 

Fz  - external  load  along  z,  lbs. 

Equation  (99)  is  a set  of  three  non-linear  simultaneous  equations  in  which 
the  variables  are  6x,  6^  and  6^.  They  are  solved  numerically  in  Appendix  A 
using  Newton-Raphson  iteration. 


M' 

x 


l [RP,  sinS. 


2(i-X.)M 

1 


. x (RcosB, 

i-1  2i  2i  , - d 2i 


n 2(1-Xi)M8l 

MV  - l [RP,  sinS, (RcosS, 

y i-1  i d Zi 


If  6,  are  current  estimates  of  the  three  displacements,  improved  estimates 
m-1 

are 


where 

V 


*■-  ■ V,  'i'.J'o 


(100) 


are  the  residues  of  Equation  (99)  evaluated  at  current  estimates  of 


The  elements  of  the  coefficient  matrix  are  the  partial  derivatives  of  the 
bearing's  reactions  with  respect  to  the  displacements  5^  and  comprise  the 
stiffness  matrix,  Equation  (13). 


The  elements  of  the  stiffness  matrix.  Equation  (13),  form  the  data  input 
for  the  rotordynamic  response  computer  program  described  in  Reference  2. 
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SECTION  III 

APPLICATION  OF  COMPUTER  PROGRAM 


The  analysis  of  Section  2.0  has  been  programmed  in  Fortran  IV  for  a 
digital  computer  and  is  suitable  for  use  on  the  CDC  6600.  A program  list- 
ing is  presented  in  Appendix  A. 


3.1  Input  Format 

Instructions  for  preparation  of  the  input  cards  are  included  as  comment 
cards  in  the  program  listing  in  Appendix  A.  To  illustrate  a typical  case 
consider  the  bearing  in  Table  1. 


TABLE  1 
SAMPLE  CASE 

DEEP  GROOVE  BALL  BEARING 


Number  of  Balls 
Ball  Diameter 
Pitch  Diameter 
Contact  Angle 
Outer-race  Curvature 
Inner-race  Curvature 
Poisson's  Ratio 
Modulus  of  Elasticity 
Ball  Density 
RPM  of  Inner  Race 
RPM  of  Outer  Race 
Radial  Load 
Axial  Load 


9 

.1875  inches 
.9252  inches 
0° 

.530 

.516 

.25 

29  x 106  lb/in. 
0.283  lb/in.3 
1000  rpm 
0 rpm 
1000  lbs. 

0 lbs. 


Figure  14  shows  a printout  of  input  for  this  case. 
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3.2  Output  Format 

Figure  15  shows  the  output  format  for  the  sample  case  of  Table  1.  The 
output  data  includes  the  internal  load  distribution  and  various  other 
9tress  and  dynamic  parameters. 

The  last  section  of  data  provides  the  complete  stiffness  matrix. 
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SECTION  IV 
DESIGN  DATA 


Three  separate  sets  of  design  charts  were  included  in  the  original  Part  IV 
report.  These  were: 

a.  Pure  Radial  Lo;  ed  Bearings  (Deep  Groove)  Contact  Angle 
0'  - 0° 

b.  Pure  Thrust  Loaded  Bearings  (Deep  Groove)  Contact  Angle 
0'  ■ 0° 

c.  Angular  Contact  Bearings  with  Axial  Preload  and  Applied 
Radial  Load  0'  ■ 25°,  15° 

As  a check  of  these  design  charts  shows  some  differences  between  the 
original  program  prediction  and  the  present  program  prediction,  the  design 
charts  have  been  recalculated  2nd  are  shown  in  Appendices  B,  C,  and  D. 

Table  2 describes  the  dimensions  and  symbols  used  for  the  deep-grooved 
ball  bearings.  Table  3 contains  information  pertaining  to  the  angular 
contact  bearings.  These  bearings  are  identical  to  those  studied  previously 
in  Reference  1. 

4.1  Radial  Stiffness  Versus  Radial  Load 

The  first  set  of  three  charts  (Appendix  B)  contains  graphs  of  radial 
stiffness  versus  radial  load.  Load  levels  are  indicated  on  thp  curves. 

The  effects  of  bearing  size  and  race  curvatures  are  illustrated  by  these 
four  charts.  In  general,  a bearing  with  curvatures  of  f^  ■ .5)0, 

" .516  is  stiffer  than  the  same  bearing  operating  with  curvatures  of 
fj  ■ fj  * .570,  for  the  same  radial  load.  Radial  stiffness  is  higher 
for  a bearing  with  a larger  bore  diameter  and/or  a greater  number  of  balls. 
Note,  for  pure  radial  load,  the  linear  relationship  between  leg  K and 
log  F . This  was  previously  illustrated  ir  Figure  6. 
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4,2  Axial  Stiffness  Versus  Thrust  Load 

The  second  set  of  eight  charts  (Appendix  C)  contains  graphs  of  axial  stiff- 
ness and  axial  deflection  versus  axial  thrust  applied  load.  Load  levels 
are  tabulated  in  Table  4 for  bearings  undergoing  a pure  thrust  load. 

A similar  observation  as  given  above  for  radially  loaded  bearings,  can  be 
made  for  the  thrust  loaded  bearing,  i.e.,  a bearing  operating  with  curva- 
tures of  f^  ■ .530,  f^  “ .516,  is  stiffer  than  the  same  bearing  operating 
with  curvatures  of  fj  » f2  ■ .570  for  the  same  axial  load.  For  all 
practical  purposes,  however,  an  average  curve  may  be  drawn  for  axial  stiff- 
ness versus  axial  load  for  all  bearing  sizes.  In  practicular,  the  bearing 
with  the  smaller  bore  and  less  balls  is  less  stiff  at  light  loads  and  more 
stiff  at  heavy  loads  as  compared  to  the  larger  bore  bearing.  There  is  an 
approximate  linear  relationship  between  log  and  log  (also  see 
Figure  8) . 


4.3  Radial  Stiffness  Versus  Radial  Load  with  Preload 

The  third  set  of  (24)  charts  (Appendix  D)  contain  graphs  of  radial  stiffness 
versus  radial  load  with  preload.  Load  levels  are  indicated  on  the  curves.  The 
effects  of  bearing  size,  race  curvatures,  initial  contact  angle,  and  axial 
preload  are  illustrated  by  these  24  charts.  For  the  same  radial  load  and 
axial  preload,  a bearing  operating  with  curvatures  of  f^  «■  .530.  ” .516. 

is  stiffer  than  the  same  bearing  operating  with  curvatures  of  * 

.570.  The  radial  stiffness  level  is  higher  for  a bearing  with  a larger 
bore  diameter  and/or  a greater  number  of  balls,  and  the  smaller  initial 
contact  angle,  (S'  “ 15°).  In  general,  the  radial  stiffness  vs.  radial  load 
curve  for  an  angular  contact  bearing  is  composed  of  three  different  behav- 
ing regions.  One  region  shows  the  stiffness  to  be  constant  with  varying 
radial  load.  (This  is  the  light  radial  load  region.)  The  middle,  or 
moderate  radial  load  region  shows  a minimum  value  for  radial  stiffness. 

The  heavily  radial  loaded  region  shows  a linear  relationship  between  log 

K and  log  F . This  third  region  is  similar  in  behavior  to  that  of  the 
XX  X 

characteristics  of  a pure  radial  loaded  deep  grooved  bearing.  The  basic 
cause  for  this  curve  having  three  separate  regions  is  due  to  the  axial 
preload.  In  region  one,  the  axial  preload  has  a great  effect  in  holding 
the  radial  stiffness  constant.  In  region  two,  where  the  applied  radial 
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TABLE  4 


AXIAL  LOADED  DEEP  GROOVE  BEARINGS 

Table  of  approximate  load  values  corresponding  to  C/P  ■ 5 and  C/P  ■ 10  load  levels. 


Bearing  Symbol Load  (lb) 


C/P  - 5 

C/P  = 10 

Al 

290 

100 

B1 

600 

250 

Cl 

1100 

450 

DI 

2250 

950 

El 

3650 

1500 

A2 

70 

30 

B2 

175 

75 

C2 

300 

100 

D2 

550 

250 

E2 

950 

350 

AA1 

380 

155 

BB1 

800 

300 

CC1 

1550 

650 

DDl 

3000 

1250 

EE1 

5050 

2100 

AA2 

95 

50 

BB2 

200 

100 

CC2 

400 

200 

DD2 

1500 

350 

EE2 

2550 

700 

UU 


load  becomes  equal  in  magnitude  to  che  axial  preload  the  radial  stiff- 
ness tends  to  decrease  with  increasing  applied  radial  load  to  a minimum 
value.  In  the  third  region,  the  axial  preload  has  little  or  no  effect, 
and  the  angular  contact  bearing  reflects  the  behavior  of  a pure  radially 
loaded  bearing  i.e.,  a linear  log  K versus  log  F relationship. 

Thus  another  point  one  is  led  to  observe  is  the  role  of  axial  preload 

magnitude  on  the  three  regions  of  a typical  stiffness  versus  load  curve. 

Three  different  preloads  are  represented  in  these  charts  and  are  tabulated 

in  Table  3.  These  preloads  are  given  the  names  selected  light  moderate, 

and  preferred  heavy.  The  effect  of  increased  preload  is  to  increase  the 

region  one  load  range  and  decrease  region  three  load  range.  Thus,  the 

ultimate  is  a constant  radial  stiffness  with  varying  radial  load  obtained 

with  an  infinite  preload.  The  increased  preload  also  has  the  effect  of 

increasing  the  level  of  stiffness  in  regions  one  and  two.  However,  it 

should  be  noted  particularly  that  the  level  of  stiffness  in  region  three. 

for  the  same  radial  load,  is  the  same  for  all  preload  values.  This,  as 

mentioned  above,  is  because  the  axial  preload  effect  is  relieved  entirely 

above  a certain  (radial  load)  (axial  preload)  ratio.  (Approximately  F / 

F = 3 for  6’  » 25°  and  F /F  =4  for  S’  «=  13°. ) 

Z X z 

In  general,  the  light  and  extra  light  deep  grooved  ball  bearings  examined 

here  will  have  a radial  stiffness  ranging  from  10'  to  2 x 10^  for  radial 

loads  of  from  10  to  2,000  lbs.  The  angular  contacl  bearings  will  have 

radial,  stiffness  values  from  2 x JO'’  to  2 x 10^  for  radial  loads  of  from 

10  to  2,000  lbs.  The  deep  grooved  ball  bearings  will  have  an  axial  stiff- 

4 6 

ness  per  bearing  of  from  2 x 10  to  4 x 10  for  thrust  loads  of  from  10  to 
4 

10  lbs.  As  in  the  case  of  the  preloaded  radial  bearing,  preloading  will 
Increase  these  values  of  axial  stiffness. 


APPENDIX  A 

COMPUTER  PROGRAM  FOR  CALCULATING  THE 
STIFFNESS  MATRIX  OF  A BALL  BEARING 
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APPENDIX  B 

BEARING  STIFFNESS  DESIGN  CHARTS 
PURE  RADIAL  LOAD 
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BEARING  STIFFNESS  DESIGN  CHARTS 
PURE  THRUST  LOAD 
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APPENDIX  E 

BEARING  STIFFNESS  DESIGN  CHARTS 
ANGULAR  CONTACT  BEARING  WITH  PRE-LOAD 
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Fig.  K-4  Radial  Stiffness  for  Angular  Contact  Bearing, 
Preload  — Selected  Light, 


Radial  Stiffness  for 
Preload  — Preferred 
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